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The perturbative three-loop result for the thermodynamic potential of QCD at finite temperature 
and chemical potential as obtained in the framework of dimensional reduction is compared with the 
exact result in the limit of large flavor number. The apparent convergence of the former as well 
as possibilities for optimization are investigated. Corresponding optimized results for full QCD are 
given for the case of two massless quark flavors. 
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I. INTRODUCTION 

Quite some progress has been made recently in extend- 
ing perturbative and nonperturbative results on the ther- 
modynamic potential of hot deconfined QCD to nonzero 
chemical potential. On the lattice, where a nonzero quark 
chemical potential gives rise to a computationally prob- 
lematic complex action, a number of new techniques have 
been devised that make it possible to study the effects of 
a not too large quark chemical potential < T, and 
first results have been produced for the equation of state 
P, 0, ■ Concurrently, the perturbative series expansion 
of the QCD pressure, which had been determined up to 
and including the last fully perturbative order In g in 
Refs. 0, IE IE 0) has been extended to non-zero chem- 
ical potential in Refs. |EJE_|^ within the framework 
of dimensional reduction [a ll^ . Dimensional reduction 
can be expected to work as long as 27rT is the domi- 
nant energy scale. In the limit of large ^/T this ceases 
to be the case when the scale gfj, appearing in the three- 
dimensional mass parameter m e becomes of the same or- 
der. In the regime T < g/j,, perturbation theory requires 
a reorganization with new phenomena such as the onset 
of non- Fermi- liquid behavior [l^ . and for the paramet- 
rically smaller temperatures T < fig~^ exp[—3Tr^ / {^/2g)] 
eventually the nonperturbative phenomenon of color su- 
perconductivity arises [l3| . 

Since at any temperature of practical or even theoret- 
ical relevance, the QCD coupling g is not much smaller 
than 1, the predictivity of the perturbative results re- 
mains uncertain. In fact, it is a well known problem 
of perturbative results at finite temperature that the 
strictly perturbative series expansion has extremely poor 
apparent convergence and large renormalization scale de- 
pendence up to temperatures of at least 10^ times the 
deconfinement temperature Tc [E • On the other hand, 
it has been found that already resummations involving 
only the lowest order contributions from the soft scale gT 
can improve the situation decisively (' jl4j and references 



therein). Within the framework of dimensional reduc- 
tion, a significant improvement can be achieved by the 
simple resummation provided by keeping effective-field- 
theory parameters unexpanded 0, Il5j . 

In particular at non-zero quark chemical potential, 
where lattice data are still somewhat preliminary, it 
would be helpful to be able to test the predictivity of the 
(resummed) perturbative results by other means. Most 
recently the limit of large flavor number Nf — > 00, with 
finite g'^g = g'^Nf/2, has been developed as such a theo- 
retical test bed. In Refs. |3| the large- iV/ pressure has 
been determined exactly for all temperatures sufficiently 
below the scale of the Landau pole A^ ~ T exp{6TT^ g^g) 
that is present in the large- A^/ theory; in Ref. [13 this 
has now been extended to non-zero chemical potential. In 
this Brief Report, we use the results of Ref. |17] obtained 
in the large- A'^y limit to investigate the apparent conver- 
gence of the perturbative results of Ref. obtained 
in dimensional reduction at non-zero chemical potential. 
We are also able to study for the first time quantita- 
tively the inevitable breakdown of dimensional reduction 
at large fj,/T in the context of the pressure. For mod- 
erate IJ,/T, we consider the possibilities for optimization 
of the perturbative results, and present correspondingly 
optimized results for Nf — 2 massless flavors in full QCD. 



II. LARGE iV 



In the limit of a large number of quark flavors Nf with 
a common chemical potential fi, the dimensional reduc- 
tion result for the interaction part of the pressure of QCD 
at finite temperature and chemical potential obtained in 



[13 reduces to {g^g = g^Nf/2) 
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where Ng is the number of gluons. The effective field 
theory parameter m|; is given by 
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the special functions [iC 
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C is the Riemann zeta function, C,'{x,y) = dxC{x,y), and 
^ the digamma function \E'(w) = r'{w)/r{w). Note that 
despite the appearance of complex quantities in ||2Jl the 
coefficient a^a is real as it has to be. 

To the order of accuracy of the result Q , the relevant 
dimensionally reduced effective theory in the large-TVj^ 
limit turns out to be noninteracting, contributing only 
the term m'^T in This gives rise to nonanalytic terms 
in g2g, namely single powers of (jcff, but no logarithmic 
terms as in full QCD. 

Just as in full QCD, however, the perturbative result 
has unusually large renormalization scale dependence for 
5eff ^ 2. But, as noted in Ref. ilfl], choosing the renor- 
malization scale p,MS such that the g^g correction in rn^ 
as given in Eq. is put to zero (FAC-m) gives good 
agreement with the exact result up to g'^g 7 at zero 
chemical potential. At large Nf, this prescription is also 
FAC with respect to the pressure itself, since it puts to 
zero the coefficient of its g^g-term. 

In the upper panel of Fig. ^ the quality of the per- 
turbative result to order g^g in the laige-Nf limit with 



fJ-MS 



Mms^ ™ displayed for the entire ^-T plane. 



The deviation of the thus optimized perturbative re- 
sult for the interaction part of the pressure from the 
exact one is shown in the form of a contour plot of 

[{Pdr-Po)~{P-Po)]/\P-Po\ = {Pbk~P)/\P-PoI 
where P is the exact result from Ref. JJj and Pq the 
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FIG. 1: Percentage errors of the perturbative result for the 
interaction part of the pressure to order g^^g in the large- A^/ 
limit for two choices of p.MS'- Fastest apparent convergence of 
P as well as m% (FAC-m), and of gl (FAC-g). The brightest 
area corresponds to an error of less than 1%, the darkest ones 
to an error of over 100%. The ratio of chemical potential to 
temperature increases from top to bottom according to <^ = 
arctan(7rr//x) so that 90° corresponds to high temperature 
and zero chemical potential, and 0° to zero temperature and 
high chemical potential. The coupling is given in terms of 



ffcff(MMs) at ^J,MS = yV2T2 + 



ideal-gas value. The fi-T plane is parametrized by 1^9 = 
arctan(7rT//z) and ^^^(/xms) at /zms = \/ tt^T^ + fj,"^. The 
accuracy of the FAC-m result at zero chemical potential 
{tp = 90°) is seen to decrease slowly with increasing chem- 
ical potential, apart from an accidental zero of the error 
when tp ~ 22°; for slightly smaller ip < 18°, i.e. T < O.l/i, 
the errors eventually start to grow rapidly, marking the 
breakdown of dimensional reduction. 

While FAC-m is particularly natural in the large- A^/ 
limit, one can consider different optimizations. At sub- 
leading orders in the large- A'^y expansion, one would also 
need the cffcctive-ficld-theory parameter 
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and one could choose to set /Ims such that the g'^g cor- 
rection in g"^ is put to zero (FAC-g). Curiously enough, 
the corresponding /2ms happens to be strictly propor- 
tional to that of the FAC-m scheme, /iMcp~^(T, /i) = 
In contrast to FAC-m, this does not 

Neverthe- 
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less, as can be seen in the lower panel of Fig.Q] the qual- 
ity of the perturbative result is rather similar to FAC-m: 
the sign of the deviation is now positive throughout the 
^-T plane so that there is no accidental zero of the er- 
ror, but the absolute magnitude of the error is otherwise 
comparable. 

Because FAC-m and FAC-g have comparable errors 
with opposite sign for tp > 18°, the optimal choice of 
renormalization scale for the 3-loop result happens to lie 
in between the two FAC scales. Incidentally, their arith- 
metic mean turns out to give a surprising accuracy with 
errors below 1% for all values of g^g considered, if ^ 67° 
or fi/T ^ 1.3. However, for larger fi/T the range of cou- 
pling with small errors shrinks quickly; for T < O.l/i the 
errors are finally no smaller than with cither FAC-m or 
FAC-g. 

Thus, except at rather small T/fi, the two FAC op- 
timizations that we have considered appear to give very 
satisfactory results in the large- iVj limit, even at coupling 
so large that the renormalization scale dependence is 
already quite strong. The renormalization scales they 
lead to are in fact rather close to each other — they differ 
just by a factor of e^^^ ~ 1.6, and they (perhaps fortu- 
itously) turn out to enclose the actual optimal choice for 
T/At> 0.1. 



III. FINITE iV 



It is of course far from guaranteed that the optimiza- 
tions FAC-m and FAC-g will work equally well in full, 
i.e. finite- A^/ QCD (in the deconfinement phase). When 
applying these prescriptions now to the case of Nf = 2, 
we shall therefore also consider the size of the renormal- 
ization scale dependence as a further criterion and inves- 
tigate whether it can be reduced by a simple resumma- 
tion consisting of keeping effective-field-theory parame- 
ters unexpanded. 

In Ref. ^3 , it has been found that at zero chemical po- 
tential such a resummation considerably reduces the scale 
dependence. Using a standard 2-loop running coupling 
one can even fix the renormalization scale by a principle 
of minimal sensitivity (PMS), i.e. by requiring that the 
derivative of P with respect to /Ims vanishes, and this 
turns out to be close to the FAC result.^ While the for- 
mal /IMS-independence of the perturbative result to order 
g5 Qj. g6 g requires only a one-loop beta function, it is 
not inconsistent to use a more accurate coupling. We 
in particular adopt the 2-loop running coupling because 
the QCD scale Aqcd of other renormalization schemes is 
frequently related to the one of the MS scheme using the 
2-loop beta function. Also, the 2-loop coupling is already 
reasonably close to the 3- or 4-loop coupling. 
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FIG. 2: Three-loop pressure in deconfined QCD with two 
massless flavors for fi = Aqcd (upper panel) and /i = 5Aqcd 
(lower panel). The bright-colored bands show the strictly 
perturbative results to order with fiMS varied about the 
FAC-m value (dash-dotted line) by a factor of 2, the darker 
bands the same when m% is left unexpanded. The dashed line 
is the FAC-g result, the full line beginning at T ~ I.SAqcd 
the PMS result. 



In Fig. 12 we have repeated the analysis of the scale 
dependence of for the perturbative Nf — 2 result of 
Ref. ^3 to order^ with a quark chemical potential of 
/i = Aqcd as well as = 5 Aqcd- The former value is 
within what is presently being achieved in lattice calcu- 
lations for Nf — 2 3] ; for /i = 5Aqcd one may expect to 
be in a quark-gluon phase for all values of the tempera- 
ture. In the plots of Fig. |21 the temperature is varied in 
both cases between 0.5 and 5 times Aqcd- 

The case fi = Aqcd turns out to be qualitatively simi- 
lar to the ^ = case studied in [l^ : the scale dependence 
is greatly reduced by keeping unexpanded, and the 
optimization using FAC is close to PMS (where the lat- 
ter exists). For T < Aqcd the scale dependence increases 
strongly so that one probably should conclude that any 
predictivity is lost. On the other hand, all optimized re- 
sults are very close to each other for T > 1.5 Aqcd and 
may be considered as rather definite predictions. 

In the case /i = 5Aqcd, the situation is similar for 
T ^ 1.5 Aqcd, but the scale dependences turn out to 
increase more rapidly as T is lowered; moreover, keep- 



^ In the large-A^^ limit PMS cannot be applied because there the 
scale dependence turns out to be monotonic. 



^ We refrain, however, from extendi ng t his analysis to the order 
\ng contribution as done in Ref. Il5l for zero chemical poten- 
tial, because it would involve an unknown function of fj,/T rather 
than a mere constant. 
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FIG. 3: The diflterence AP = P(T, /x) - P(r, 0) divided by 
using the unexpanded three-loop result from dimensional 
reduction of Ref. 10] for [ijT = 0.2, . . . , 1.0 (bottom to top). 
Shaded areas correspond to a variation of /^ms around the 
FAC-m choice by a factor of 2; dashed and dash-dotted lines 
correspond to the FAC-g and FAC-m results, respectively. 
Also included are the recent lattice data of Ref. Q (not yet 
continuum-extrapolated!) assuming T^^" = 0.49 Aqcd [iq- 

ing to|; unexpanded does no longer help to improve this. 
Yet, the two FAC results are fairly close to each other 
for T ^ O.SAqcd, but show a run-away behavior for 
T slightly smaller than O.SAqcd (not displayed). This 
could indicate a beginning breakdown of dimensional re- 
duction, or at least of its predictivity. On the other 
hand, the FAC results appear to remain predictive for 
T > 1.5Aqcd even when [i ~ SAqcd- At much higher /x, 
however, T needs to be correspondingly higher in order 
that reasonable FAC results are obtained, but this is pri- 
marily due to the run-away behavior of the FAC scales 
and not so much due to the behavior of the pressure itself. 

Another popular choice of renormalization scales is the 
BLM prescription 19], which amounts to simply keeping 
the FAC scales as obtained in the iV/ — > co limit. These 



scales show no run-away behaviour, but are significantly 
smaller than those of FAC-g and FAC-m at finite A^/, 
with results roughly at the lower boundaries of the bands 
of Fig. 13 i.e. also comparatively far from PMS. 

We also note in passing that just as in the N f — > oo 
limit the two FAC scales happen to be simply propor- 
tional to each other when Nf — 2: Mms^'^/^ms^'™ = 
g5/29 ^ 1.188^ while for N f ^ 2 this ratio has nontrivial 
/i/T-dependence. 

In Fig. Owe display the optimized results for the differ- 
ence AP — P{T, fi) — P{T, 0) for various fi/T correspond- 
ing to the recent lattice results given in Fig. 6 of Ref. 
assuming Tq = T^=° = 0.49 Aqcd "T^- At T/Tq = 2, the 
highest value considered in ^ , our FAC-g and FAC-m re- 
sults exceed the not-yet-continuum-extrapolated lattice 
data consistently by 10.5% and 9%, respectively. This 
is in fact roughly the expected discretization error |2Clj |. 
When normalized to the free value APq instead of T"', the 
results would be essentially /^-independent and thus also 
very similar to the TV/ = 2 + 1 lattice results of Ref. Ql as 
well as the quasiparticle model results of Refs. 0, |23| ■ 

To summarize, we have found that the requirement of 
fastest apparent convergence of effective-field-theory pa- 
rameters works remarkably well when the perturbative 
3-loop result for the pressure of QCD at finite tempera- 
ture and chemical potential is compared with the exactly 
solvable large- A''/ limit, except for the region T < 0.1/x. 
The perturbative results for finite Nf turn out to agree 
reasonably well with existing lattice data for deconfined 
QCD with non-zero fj,, and also the otherwise strong 
renormalization scale dependence is brought under con- 
trol, when effective-field-theoryparameters are kept un- 
expanded as proposed in Refs. 
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